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Signal  detection  problems  are  treated  in  the  context  of  several  Gaussian 
and  non-Gaussian  stochastic  processes  with  one  or  more  nuisance  parameters. 
The  techniques  all  entail  a  basic  data  transformation  and  the  use  of  a 
Kolmogorov-type  statistic  on  one  or  more  components  of  the  transformed  data. 
Extraneous  statistical  noise  and  randomized  rank  statistics  are  employed. 
Numerical  examples  are  given.  ^ 
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1-  Intferitoation  JEA  Swmry. 


One  decision  problea  encountered  in  signal  detection  is  the  situa¬ 


tion  in  which  the  law  (L)  of  the  process  in  the  pure  noise  (PN) 


ease  is  known  precisely  (i.e.,  L  .*  1*),  and  in  the  noise  plus  signal 
»(N  ♦  S>  ease  L  f  L*.  ,  Me  assune  the  data  consists  pf  N  observations 
Xj,  ,  XN  of  the  process.  For  this  situation  standard  goodness -of- 
fit  (G-O-F)  detectors,  sudl  as  those  based  on  the  Kolaogorov-Bnixnov 
and  Craner-von  (Uses,  statistics  can  be  used. 

One  considers  here  extensions  of  theabove  problea  to  the  case  diere 
the  pure  noise  law  is  known  up  to  sons  nuisance  paraaaters.  The  detectors 
are  all  bated  on  the  anginal  statistical  noise  CM-S-N)  whidi  is  dis¬ 


tributed  independently  of  the  ainiaal  sufficient  statistics  (M-S-S) ' 
and  togetlMr  with  the  M-S-S  constitute  a  one-to-one  (a.e.)  sapping  of 
the. date.  Ibis  trws  formation  .«  called  the  basic  data  trsnsfbmation 
(IDT) .  The  M-S-N  and  its  u  tion  fbr  several  copaon  processes 


axe  presented  in  Section  2. 


Lilliefors  (1907,  1909)  has  developed  procedures  for  the  noraal  and 
exponential  oases  when  nuisance  psrsaeters  ere  present  by  using  the 
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maximum  likelihood  estimate  of  the  distribution  function.  For  the 
Kolmogorov- Smirnov  statistic  Srinivasan  (1970)  has  developed  similar 
procedures  based  on  the  Rao- Blackwell  estimate  of  distribution  function. 
Their  results  are  presented  in  Section  3  along  with  extensions  to  the 
uniform  distribution  done  by  Choi  (1980).  Also,  results  for  the  normal 
case  when  the  mean  is  known  and  the  only  nuisance  parameter  is  the 
variance,  are  presented. 

Section  4  illustrates  how  several  other  common  detection  procedures 
can  be  applied  to  the  M-S-N  when  there  are  nuisance  parameters.  The 
types  of  procedures  applied  therein  include  (a)  G-O-F  (e.g.,  Cramer- 
von-Mises),  (b)  classical  (e.g.,  F-statistie) ,  and  (c)  nonparametric 
(e.g.,  sign-rank  statistics). 

The  use  of  extraneous  statistical  noise  (E-S-N)  in  detection 
procedures  is  introduced  in  Section  5.  E-S-N  is  a  simulated  sample 
with  known  law  LQ  belonging  to  the  same  family  as  l*.  By  applying 
the  inverse  of  the  BUT  one  obtains,  under  PN,  a  sample  with  law  LQ. 
One  can  then  apply  standard  detectors  for  this  known  distribution.  The 
section  concludes  with  a  table  of  the  inverse  of  the  BUT  for  several 
stochastic  process  models  of  interest. 

Finally,  in  Section  6  all  of  the  above  mentioned  detection  techni¬ 
ques  are  applied  to  computer  generated  data  sets  of  five  different 
processes.  In  addition,  the  techniques  are  applied  to  two  real  data 
sets  consisting  of  interarrival  times  of  epileptic  seizures  for  two 
epileptic  females. 


\  2.  Statistical  Preliminaries:  Maximal  Statistical  Noise  (M-S-W). 

Many  decision  procedures  are  based  on  the  minimal  sufficient 

statistic  (M-S-S).  This  is  particularly  so  when  the  inference  is  centered 

on  the  value  of  a  parameter.  However,  in  this  paper  one  is  primarily 

concerned  with  inference  in  the  presence  of  nuisance  paraaeters.  Hence, 

the  need  to  utilize  s one thing  other  than  the  M-S-S.  This  other  quantity 

will  be  called  the  naxinal  statistical  noise  (M-S-W). 

Definition  2.1  -  Let  S(X)  be  the  M-S-S  for  data  X  *  (Xj ,  ...  t  X^) , 

where  the  law,  L,  of  the  tiae  series  is  an  eleaent  of  the  faaily  ft' . 

Further,  let  S(X)  »  [N(X),  S(X)]  be  a  traits foraation  which  is  1-1 

a.e.  and  where  N(X)  and  S(X)  are  statistically  independent.  Then 
*  ^  ^ 

(a)  W(X)  is  called  M-S-W,  i.e.,  aaxiaal  statistical  noise;  and 

‘Vi 

(b)  6(X)  is  called  the  BDT,  the  basic  data  trans foraation. 

Bxaaple  2.1  -  (Homogeneous  Poisson  Process,  HPP) 

Let  [N(t):  t  _>  0]  be  a  HPP,  i.e.,  hoaogeneous  Poisson  Process.  Let 

(t  )  and  (W  )  be  the  associated  interarrival  tiae  and  waiting  tiae 
n  n 

series,  respectively.  If  one  views  the  data  as  X  ■  (Wj,  ....  *n),  then 
one  has 

(a)  M-S-S:  S(X)  •  WM 

H  N  If  If 

(b)  M-S-N  N(X)  -  (Y1#  ...,  YN1)  -  (^  ,  ^  ,  ....  X,  and 

^  N  N  N 

(c)  BDT:  «({)  -  [«(X),  S(X)J. 

Bxaaple  2.2  -  (Wiener- Lew  Process,  WLP) 

4  Let  W(t)  ■  u(t)  ♦  0W*(t),  where  [W*(t):  t  >  0]  is  a  Wiener- Levy' 


(WL)  process  satisfying  (i)  I{W*(t)>2  0;  (ii)  Cov  (lf*(s).  »*(*)>  - 
■in  (s,t);  (iii)  it  is  Gaussian.  Let  Zf  -  W(rA),  i.e.,  eat  sables 
st  tiass  b,  2b,  3b,  ,  and  X  ■  (Xj,  ...»  X^)  where  Xy  ■  Zy  -  Z^ 

and  ZQ  >0. 

Casa  I.  u(t)  s  gt 

Hare,  Xj,  ....  X#  are  i.i.d  M(8A,  a2A); 

+  II  k 

M-S-S:  S(X)  •  (I,  S%)  where  X-H"1  J  X^,  S2  »  IT1  |  CX^-IF)2 

X.  -  J  K.  -  X 

M-S-M:  M(X)  •  C-i* - ,  ...»  -V~  );  and 

*  *X  °X 

BUT:  «(*>.  ptaf),  Sft)] 

Casa  II.  ltftl  »  0 

Here,  Xj,  ....  X„  ara  i.i.d.  H(0,  o2A)i 

M-3-S:  SQC) 

M-J4I:  *(X)  •  (Bj,  ....  Vj,  ...»  VM);  what* 

tT  •  «(Xy)  with  c(U)  •  1  if  U  >  0  and  -  0 


if  U  <  0; 


[One  should  note  here  that  if  "o"  is  a  nondegenerate  randoa  variable, 
then  {h(t) }  is  non-Gauss ian,  but  many  of  the  detections  techniques 
for  Gaussian  models  still  apply.] 


Example  2.3  -  (Uniform  Renewal  Process,  URP) 

Let  Uj,  ...,  Ujj  be  i.i.d.  U(0,8)  be  the  interarrival  times  of  a 
(URP)  point  process.  Then  one  has 


M-S-S:  S(U)  •  U(N); 

a* 


M-S-W:  «(U)  •  (R,  V),  where 


$  •  [RCUj),  ....  R(Un)]  and  V  -  [■ 


ua 


U 


mn 


LL5I £1  1. 

ufiir 


BUT:  Sty  -  [W(X),  S(X)] 


Example  2.4  -  (Non-homogeneous  POisson  Process,  NHPPj) 

Let  [N*(t):  t  _>  0]  be  a  FMPP  whose  mean  function  is  known  up  to  a 


multiplicative  (nuisance)  parameter  a,  i.e.,  u*(t)  ■  a  yQ(t),  and 


let  . W*  •  (W*,  ...»  Wj*)  be  the  waiting  time  data. 

-1. 


Now  form  M(t)  ■  N*[uQ  (t)  3 .  Then  [M(t)]  is  a  HPP  with  mean  function 
U(t)  *  at  and  waiting  times  W  »  (Wj,  ...,  W^)  where  *r  »  U0(W*). 

Hence,  one  has 


M-S-S:  S(W)  -  VqOIJ); 


I  /con^lon  For 


M-S-N: 


N(W)  •  (V.,  ...»  V„  .)  with  V.  - 


Uq  TO 


-V  '1* . N-l'  'r  '  U^T^T 


g::  \&i 

i.y', 

;o  r.'.;  cod 


■i.  I?  lection. 


’  •  ■*  n:*!  /  1 


BUT: 


•  [N(W),  sqp). 


The  statistics  in  the  examples  above  involve  the  fol lowing  distributions 
Definition  2.2  -  (Basic  Distributions) 

(1)  If  xj»  •••»  \  are  i‘i‘d  F(.) »  continuous,  one  says 

PtfXj).  ....  RfX^]  *  R  -  V(k). 

(2)  If  Oj,  ....  ere  i.i.d.  U(0,1)  then 

[U(l),  ....  U(k)J  -v  U-O-S0O . 

(3)  Let  X1#  ....  Xfc  be  i.i.d.  N(y.o2),  then 

X-J  X.  -  X 

iMLS-N(k),  where 

5X  5X 

|*r  | 

(4)  If  P(V  <  x)  ■  xk,  0  <  x  <  1,  then  one  says  V  *  PW(k). 

These  processes  and  statistics  are  suflitrixed  in  Table  2.1  below 
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0  <  2  <  Y(n) , 


(viii) 


G*(9;i) 

G**(0;z) 


■  U(0,Y(n)) 


Y(n)  ’ 

Z  >  Y(n). 


Then, 

(1)  Dj  *  sup  |Fy  ^(z)  -  G(z)|  'v  K  -  S(k),  i.e.,  has  the  K-S  distribution 


for  a  random  sample  of  size  k  (See  Bimbaua,  19S2). 


(2)  D2  *  sup  |F*k)  -  6(A;z)|  %  LLEX(k)  if  H*  Exn  (A)  (See  Lilliefors,  1970). 

(3)  Dj  «  sup  | -  G(A;z)|  %  SREX(k)  if  G  »  Exn  (A)  (See  Srinivasan,  1970). 

z  7 

(4)  D4  -  sup  |F^k)(z)  -  G(p;a2,z) |  ^  LLGA(k)  if  G  -  N(u,o2)  (See  Lilliefors, 
(1969) . 

(5)  Ds  -  sup  |F^k)(z)  -  G(jj,o2;z)|  -v  SRGA(k)  if  G  *  N(y,o2)  (See 

z  7 

Srinivasan,  1970). 


(6)  D6  .  sup  |F^k)(z)  -  G*  (0;z)|  -v  CHL(k) ,  if  G  «  U(O,0)  (See  Choi,  1980 

z 

and  Table  A.l). 

(7)  D?  »  sup  |F^k)(z)  -  G**  (0;z)|  ^  CHS(k),  if  G.  U(O,0)  (See  Choi,  1980 

z 

in  which  it  is  proved  that  D?  ^  (k  ^  *)K  -  S(k  -  1).) 

The  last  six  statistics  are  useful  when  the  form  of  G  is  known, 
but  the  value  of  the  (nuisance)  parameter  is  unknown.  These  are  based  on 
the  works  of  Lilliefors  (1967,  1969),  Srinivasan  (1970),  and  Choi  (1980). 


who  give  the  relevant  tables.  S.  M.  Lee  computed  an  improved  version  of 
Choi's  original  table.  It  is  Table  A.l  of  the  appendix. 

The  utilization  of  these  statistics  for  the  detection  models  of  this 
paper  can  be  summarized  in  the  table  below,  where  usage  of  randomized 
noise  is  indicated.  This  concent  will  be  defined  in  the  next  section. 


Table  3.1.  GOF  Statistics 


Stochastic 

GOF 

Use  of 

Process 

Statistics 

Randomized  Noise 

URP 

CHL,  CHS,  K-S 

Yes,  with  K-S 

WLPL 

LLGA,  SRGA 

Yes,  with  K-S 

U(t)  =  lit 

NLP 

LILE,  SRLE 

Yes,  with  K-S 

■C 
/— \ 

rt 

w 

III 

O 

(See  discussion  below) 

HPP 

HEX,  SREX,  K-S 

Yes,  with  K-S 

u(t)  =  At 

NHPP 

LLEX,  SREX,  K-S 

Yes,  with  K-S 

u(t)  *<*w0(t) 

The  GOF  statistic 

for  the  WLP  with  y(t) 

=  0  has  not  been  studied 

by  the  above  named  authors.  However,  one  can  derive  the  appropriate 
statistic  in  a  manner  parallel  to  those  of  Lilliefors,  Srinivasan  and  Choi. 


The  Lilliefors  idea  is  to  replace  the  unknown  parameter  by  its  MLE 

(maximum  likelihood  estimate)  in  the  distribution  function.  Srinivasan's 

approach  is  quite  different  but  s one times  asymptotically  equivalent.  His 

fk! 

idea  is  to  replace  the  empirical  distribution  function  F^  '(•)  by  its 
Tao-Blackwell  estimate.  (See  e.g.,  Bickel  and  Doksm,  1977). 


Definition  3.2. 

A 

Let  (1)  G(o;z) 

S' 

G(o;z)  * 


-  ♦(«£  [I  y2]"1/2);  and 
1 

k  k 

1  -  1/2  F(k-l,l;  (l  X2  -  z2) (N-l)”1z~2)  for  0<z<(jx2)17 

1  1 

1  if  z  >  (][  X2)1/2 
1  r 


% 

with  G(o;z)  ♦  G(a;-z)  ■  1  for  all  z. 

Here  F(m,r; •)  is  the  cdf  of  Fisher's  F-statistic  with  m  and  r  degrees 
of  freedom. 

Further  define 


(3)  Da  -  sup  IF^C*)  -  8(o;z)|  *  LILBfk)  if 

8  z  7 

G  ■  N(0,o2);  and 

% 

(4)  Dq  -  sup  iF^fz)  -  G(o;z)  |  SRLE(k)  if 

*  z  7 

G  -  N(0,O2) 


Tables  for  these  two  statistics  are  given  in  the  appendix. 

It  can  be  proved  that  the  K-S  statistic  is  a  function  of  the  M-S-S 


when  the  parameters  are  completely  specified.  The  other  eight  statis¬ 
tics  above  satisfy  a  different  property. 


Proposition  3.1  -  The  GOF  statistics  in  (2)  -  (7)  of  Def .  3. 1 
and  (3)  -  (4)  of  Def.  3.2  are  functions  of  the  data  solely  through 
their  respective  M-S-N’ s  (The  proof  will  only  be  given  in  2  cases.) 


Proof  for  Gfu.o  :z) 


K  _  It 

sup  If  I  e  (z  -  Yr)  -  ♦(■Zg~-Y)  |  -  sup  l  e(w  - 
z  1  y  w  l 

Y  -  Y 


-  •(*)!, 


where  Vy  » 


Proof  for  GfX;z) 


*up  _  l  e(z  -  Y  )  -  1  ♦  (1  -  )k-1| 

0<z<kY  K  1  T  kY 


-  sup  |l  -  (1  -  V)*’1  -l  f  CU  -  V  *  V  ) | 
0<U<1  1  r  r  i 

h 

where  V  ■  —  -  ■  ■  -  . 
r  kY 


4.  Generalized  Randomized  Rank  Procedures 

In  order  to  avoid  certain  problems  with  distributions  of  statistics 


of  interest,  Durbin  (1961),  Bell  and  Doksua  (1965)  and  others  introduced 
extraneous  noise  into  the  decision  procedures. 


Example  4.1  -  Let  Z  »  (Z.,  . Z  )  be  i.i.d.  with  continuous 

i  n 

cdf  F(.).  Let  S  *  (Sj.  ...»  Cjj)  be  i.i.d.  *  and  independent  of  the 

data  Z.  Now  define  S'*  (S',  ...»  ?')  where  S'  is  a  permutation  of  S 
^  n.  l  w  ||  *  N  <x, 

satisfying  S£  -  5(1^)),  i.e.,  f  e(S£  "  “  I  e(zk  ‘  V*  Bell 

and  Doksua  (1965)  then  prove  that  S"  ^  S«  Consequently, 

i  |  «,«j»  -  rhr  1  «RPj»  *  ;  ♦  rh?- 

In  order  to  formalize  and  generalize  this  procedure,  one  needs  the 
following  definitions. 

Definition  4.1  -  Let  (rj,  ...,  rN)  be  some  permutation  of  the 

integers  1,  ...,  N,  and  let  v  ■  (v_,  v„)  be  an  N- vector.  t* 

<V  X  N  H 

is  called  the  randomized  rank  transformation  (RRT) ,.  When 

TjJ(r^,  ...,  Tjj,  Vj,  . ..,  v^)  ■  {v(tj) ,  ...,  v(Tjj)},  *  v* ,  i.e., 

N 

is  a  permutation  of  v  such  that  l  e(v(r.)  -  v  )  ■  r..  for  1  <  k  <  N. 

*  s*l  *  *  -  - 

Besides  the  RRT  defined  above,  the  procedure  in  Ex.  4.1  involves  an 

interchange  of  MSS's.  Por  that  example  the  MSS's  are  the  order  statistics. 

Definition  4.2  -  Let  X  ■  (Xj,  ...»  X^)  and  S  *  (£j#  •••»  Cjj) 
be  independent  initial  segments  of  tine  series  with  laws  L  and  L*, 
respectively,  both  in  the  family  fl'  of  distributions.  Let  6(X)  ■  (N(X), 
S(X)>  be  the  BDT  and  £'  ■  <5  {N(X),  S(S)).  Then,  if  £  is  the  original 

data. 


(1)  £  is  called  extraneous  statistical  noise*  ESW,  and 

<v 

(2)  £'  is  called  randoaized  statistical  noise,  RSW. 

(Note:  L*  is  chosen  to  be  convenient  and  tractable.) 


The  principal  result  in  this  direction  follows  froa 


Leaaa  4.1  -  Let  Y  ■  (Yj,  ...»  Y^)  have  law  Lett';  and 

(a) 

d  .i  d 

(b)  S.  *  S(Y) .  Then  6  l(N.,  S.)  *  Y. 


Adapting  the  work  of  Durbin  (1961) ,  one  has  the  exaaple  below. 


Exaaple  4.2  -  (WLP  with  u(t)  ■  ut) 

Let  X  -  (Xj,  ...»  Xjj)  be  as  in  Ex.  2.2  and  let  £»  (q,  ...»  q)  be 
i.i.d  4  and  independent  of  X.  On  has  then  that 


X  -  X  X^.  -  X 

(1)  «(X)  -  (— ^ -  .  ....  - ,  X,  Sx);  and 


x  -  X 

(2)  d'1  (A 


C  9  •  •  •  9 

5X 


x-x 

—  ,  T,  S^)  -  £'  •  (q,  ...,  q).  idiere 


X  -  X 

q  -  X  ♦  ( - § - ) .  Further, 


d 

(3)  £•  *  £,  and 


sup  |f£°(z)  -  *(z)|  Mt.SCN). 
z  4 


An  exaaple  for  non-hoaogeneous  Poisson  processes  is  as  follows : 


As  in  Ex.  2.4,  let  W  ■  (Wj,  . ..,  W^)  be  the  waiting  tines.  Then 


2  *?  .  W, 

son  -  oe  *  iu;  noo  •  t-i — 1 .  .... 

^1  Iw  Ft  M 

WN  M 


*H-1  *  *H-1 

*2*»n 


■).  Now,  choose 


« ■  «i .  Cj|)  to  be  the  waiting  tine  of  a  HPP  with  X  •  2.5,  and 


independent  of  W. 


W2  ♦  W 
r 


Further,  let  €'  -  ONflO.  S(C»,  i.e.,  £•  -  (N*>(-4 - £) 

•V  \  *v  *  k*  x  if 

N  "N 


Then  V  "  C;  and  C;  «,•  -  C;)"1  (2-j-1-)  *  PC2"*  ™  ~  2«0 


One  conpletes  this  section  with  a  table  of  6  for  the  relevant 


stochastic  processes. 


«(x)  (See  Table  2.1) 


6.  Numerical  Examples 

To  illustrate  the  use  of  the  detectors  listed  in  the  previous  sections 
data  sets  from  five  stochastic  processes  were  computer  generated.  The 
five  processes  are: 

(1)  HPP  with  rate  paraaeter  X  »  5. 

(2)  NHPP  with  aean  function  u(t)  ■  oyQ(t) 

2 

where  vQ(t)  ■  t  ♦  t  and  a  ■  5. 

(3)  WLP0,  the  Wiener- Levy'  process  with  aean  function 
p(t)  =  0  and  o2  -  5,  A  »  2. 

(4)  WLPL.  the  Wiener-Levy'  process  with  aean  function 

2 

y(t)  ■  6t  where  B  *  5,  a  >5  and  A*  2. 

(5)  URP  with  6  -  S. 


In  each  case  the  number  of  observations  is  20.  The  waiting  times 
for  the  five  data  sets  are  listed  and  plotted  in  Figures  6.1  to  6.5. 

Each  data  set  was  used  with  its  respective  pure  noise  type  of 
detectors  as  well  as  with  one  alternative  group  of  detectors.  The  false 
alarm  rate  (FAR)  was  set  at  .01  for  all  the  tests.  The  results  are 
shown  in  Tables  6.1  to  6.5. 

In  addition  to  the  simulated  data,  two  real  data  sets  consisting  of 
interarrival  tines  of  epileptic  seizures  was  obtained  from  Choi  (1980). 
The  first  data  set,  denoted  Epilepsy  I,  contains  30  observations  on 
an  el0it  year  old  epileptic  female  recorded  froa  9:00  AM  to  9:00  PM. 
The  other  data  set.  Epilepsy  II,  contains  20  observations  on  a  twelve 


year  old  epileptic  fenale  recorded  from  7:02  AM  to  7:02  PM.  Hie  waiting 
tines  are  listed  and  plotted  in  Figures  6.1  to  6.7  of  the  appendix. 

All  S  sets  of  detectors  were  applied  to  the  Epilepsy  I  data  and 
the  HPP  and  NHPP  detectors  were  applied  to  the  Epilepsy  II  data  with 
FAR  ■  .01.  The  results  are  shown  in  Tables  6.1  to  6.5  of  the  appendix. 
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Table  A.1  Critical  Values  of  D6,CHL(k)  (See  Section  3) 
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(Table  computed  by  S.  H.  Lee) 


Table  A. 2  Critical  Values  of  D«.LILE(k)  (See  Section  3) 


0.8369 

0.7975 

0.7335 

0.6603 

0.6071 

0.5663 

0.5337 

0.5048 


0.4659 

0.4438 

0.4278 

0.4095 

0.3973 

0.3872 

0.3769 

0.3690 

0.3567 

0.3479 

0.3369 

0.3344 


0.3189 

0.3172 

0.3082 

0.3054 

0.3016 

0.2903 

0.2897 

0.2843 


0.8202 

0.7217 

0.6210 

0.5665 

0.5140 

0.4766 

0.4497 

0.4250 


4005 
0.3866 
0.3709 
0.3577 
0.3411 
0.3310 
0.3215 
0.3130 
0.3024 
0.2966 
0.2896 
0.2820 
0.2781 


0.2704 

0.2627 

0.2558 

0.2551 

0.2500 

0.2444 

0.2416 

0.2356 


0.7968 

0.6566 

0.5650 

0.5074 

0.4639 

0.4288 

0.4048 

0.3823 


0 

0.3458 

0.3318 

0.3198' 

0.3067 

0.2974 

0.2875 

0.2819 

0.2712 

0.2655 

0.2598 

0.2521 

0.2473 


0.2415 

0.2361 

0.2306 

0.2266 

0.2244 

0.2189 

0.2169 

0.2121 


0.7697 

0.5992 

0.5289 

0.4695 

0.4289 

0.3994 

0.3739 

0.3530 


0.3201 

0.3069 

0.2964 

0.2838 

0.2753 

0.2674 

0.2599 

0.2511 

0.2455 

0.2408 

0.2333 

0.2284 


0.2233 

0.2183 

0.2139 

0.2095 

0.2076 

0.2031 

0.2001 

0.1975 


0.7393 

0.5507 

0.4985 

0.4403 

0.4002 

0.3752 

0.3511 

0.3318 


0.2999 

0.2886 

0.2774 

0.2656 

0.2589 

0.2513 

0.2439 

0.2356 

0.2311 

0.2260 

0.2193 

0.2141 


0.2097 

0.2048 

0.2015 

0.1968 

0.1943 

0.1906 

0.1879 

0.1853 


0.7062 

0.5294 

0.4701 

0.4182 

0.3790 

0.3534 

0.3322 

0.3137 


0.2836 

0.2725 

0.2622 

0.2513 

0.2450 

0.2374 

0.2301 

0.2226 

0.2183 

0.2137 

0.2075 

0.2023 


0.1986 

0.1938 

0.1906 

0.1863 

0.1838 

0.1803 

0.1782 

0.1748 


0.6700 
0.5107 
0.4426 
0.3969 
0.3615 
0.3355 
0.3159 
0.2981 
.2825 
0.2695 
0.2586 
0.2491 
0.2387 
0.2328 
0.2261 
0.2186 
0.2115 
0.2078 
0.2029 
0.1973 
0.1927 


0.1888 

0.1847 

0.1809 

0.1774 

0.1750 

0.1718 

0.1689 

0.1662 


(Table  computed  by  S.  N.  Lee:  20,000  repetitions) 


Table  A. 3  Critical  Values  of  D0,  SRLE(k)  (See  Section  3) 
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